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Abstract 

Momentum transfer between matter and electromagnetic field is analyzed. The related equations 
of motion and conservation laws are derived using relativistic formalism. Their correspondence to 
various, at first sight self-contradicting, experimental data (the so called Abraham-Minkowski con- 
troversy) is demonstrated. A new, Casimir like, quantum phenomenon is predicted: contribution 
of vacuum fiuctuations to the motion of dielectric liquids in crossed electric and magnetic fields. 
Velocities about 50nm/s can be expected due to the contribution of high frequency vacuum modes. 
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Electromagnetic radiation possesses energy, linear and angular momenta like any ordi- 
nary material object. As a consequence the light can exert mechanical forces on matter 
during the interaction. The forces associated with exchange of the linear and the angular 
momenta of light with material bodies were successfully measured already hundred years 
ago . However, the fundamental question of the momentum associated with a photon 

in an optically dense medium is still under discussion f3], despite that it was formulated in 
the beginning of the previous century. This problem arises from the discrepancy between 
Minkowski's!^ Gm = / d^xD x B and Abraham'sj^ Ga = 4~: / d?xE x H results, where 

H, D and B are electric and magnetic fields and inductances correspondingly. Their dif- 
ference is significant: while Minkowski's moment is directly proportional to the refractive 
index of the medium, the moment of Abraham possess inverse proportionality. 

Minkowski momentum is considered by many as unacceptable, although it was shown by 
Jones, Richards and Leslie ( JRL) that the recoil force of the light on a mirror immersed in 
liquid is directly proportional to the refractive index of the liquid j?!. This experiment was 
conducted twice with a 20 years periodlsf. However most of the theoretical works are in 
favor of Abraham's expression (for a review, see refs. |9||inl|). 

Abraham's momentum can be derived from the Poynting energy flow vector S = -^E x H 
under assumption that all energy is purely electromagnetic and relates to the mass through 
the U = mc^ relativistic formula. It corresponds to the relativistic requirement for direct 
proportionalitvpf the energy and momentum flows (the symmetry of the electromagnetic 
stress tensor) Abraham's result is also supported by statistical physics approaches 
However, as far as we know, there are no experimental data that demonstrate the inverse 
dependence of the radiation pressure on the refractive index. This, at least, allows to 
conclude that the measured momentum is not purely electromagnetic. 

The effective momentum of a photon in a dielectric medium consists of electromagnetic 
momentum and associated motion or even radiation of matter lo^. However separate identi- 
flcation of different parts proved itself to be non-trivial and sometimes led to contradictions 
with experimental data^]. Blount [3] and NelsonQ] developed a Lagrangian formalism of 
the problem, using heuristic and microscopic averaging approaches correspondingly. They 
signiflcantly clarifled the picture by associating Abraham's expression with electromagnetic 
momentum and Minkowski's momentum with phonon like pseudomomentum. Still several 
questions remained open, especially a small discrepancy of Abraham's momentum with the 
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expression derived in refs. H0- 

In this Letter the related Lagrangian and corresponding equations of motion are derived 
using relativistic formahsm. In the case of liquid dielectric, interaction of the electromagnetic 
field with matter causes motion of the latter. Thus while Abraham's expression is indeed the 
momentum of the field, the measured momentum also includes the matter contribution and 
its value coincides with Minkowski's result. Afterwards the possible vacuum contributions 
to the motion of the matter are considered. Each electromagnetic mode possesses finite 
momentum, even in its ground state Thus modification of the modes by matter can alter 
the momentum of the vacuum. The latter generally vanishes due to counter propagating 
modes that cancel each other's contribution. This situation can be different however in 
materials that are temporally and spatially asymmetric. 

The electromagnetic field in an optically dense medium is described by the Maxwell 
equations: 

^^^4f •^" = «-^^^ = 4f ■^^ = «- « 

The free electric E and magnetic B fields exist both inside and outside the matter. The 
matter response to the radiation is taken into account through derived fields D and H, 
related to E and B by the dispersion relations. In the case of linear, non- dispersive dielectric 
medium they are given by: 

D = eE,H = -B, (2) 

where e and /i are the dielectric and magnetic constants of the matter correspondingly. 

The Lagrangian approach simplifies the investigation of light-matter interaction. The 
simplicity comes from the universality that allows to use the same approach to the system 
consisting of significantly different subsystems and the natural to the Lagrangian approach 
definition of the conservation laws. The Lagrangian which is equivalent to eqs. ((T)) and ^ 



^Field 



d'-^ i"^E^ -l.B-]- (3) 



47r V2 2^ 

The first pair of the Maxwell equations corresponds to the equations of motion: 

d dLpield _ dLpield d dLpield _ dLpield 

dt~dW ~ ~dA~'dt~dW' ~ ' 

at at 



(4) 



while the second pair of (0) follows from definitions of the vector A and the scalar $ 
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potentials: 

E = —^-V^,B = VxA. (5) 
c ot 

The motion of the matter and especially its influence on the electromagnetic fleld must 
be taken into account in a combined matter- fleld Lagrangian Lmf- The linear dispersion 
relations change in moving media to: 

D = eE + ^^~^ vxB,B = fiH+ ^^~^ Exv, (6) 

C/i 



where flrst order v/c terms were taken into account [l9| 20[. They follow from relativistic 



requirements and can be derived using the flrst order Lorentz transformations: 

E ^ E + -V X B,B ^ B + -E X V, (7) 
c c 



D ^ D + -V X H,H ^ H +-D XV, (8) 

c c 

relative to (j2I)- These transformations can be applied directly to the Lagrangian Q. Sub- 
stituting (jZj) into (jni), keeping the flrst order v/c terms and adding pv'^/2 one obtains: 



Since the liquid is assumed to be incompressible, it is described by its density p and local 
velocity v only. The equations of motion of are identical to the Maxwell equations 
with dispersion relations ©• The last term of Q can be rewritten in an interaction Aj 
form, where the current j is given hj V x [E x v) {ep — 1) /fic. The latter, at least for the 
non-magnetic /i = 1 case, can be obtained by microscopic averaging procedures jl^ |ll | . 

The Lagrangian Q is explicitly independent of the space coordinates x, due to the 
homogeneity of space. Thus according to Noether theorem, the momentum [isl : 

dL dL dAi^ 
Ml 

at 

is conserved. Substituting Q into (fTUj) one obtains: 



G = J d^x \^pv + — \^-E X B- ^ E xBjj= j d^x ^pv + —E x Hj . (11) 

The corresponding angular momentum I = x x G becomes: 

1= d^x (x X pv + X E X Hj . (12) 



Therefore the conserved hnear and angular ()12|) momenta consist of the matter and 
the Abraham's field terms. The correspondence between the conserved and the measured 
momenta follows from the analysis of the Lorentz force acting on material objects 
The pv term can be obtained from the liquid's equation of motion: 

d_dL^dL 

dtdv dR' ^ ' 

where R is the matter coordinate. Far from the boundaries, OL/dR can be neglected, leading 
to: 

^ ^_t^E X B. (14) 
A-npc 

This expression corresponds qualitatively to the pseudomomentum of ref . . Substituting 
(fT^ into (fTT|) . one obtains G = D x B/Anc, which was observed in the JRL experiments. 

By analogy, the measured angular momentum is I = x x [D x B). It can be separated 
into "spin" D x A and "orbital" ^ Dj {x x V) Aj parts. The latter does not vanish only in 
the case of beams with specific wavefront distortions, e.g. Laguerre- Gaussian beams In 
contrast to the linear momentum, the spin part of the nearly plane wave is independent of 
dielectric properties of the medium. This was verified experimentally for microwave radiation 
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The dielectric constant dependent angular momentum was observed inside a cylindrical 
capacitor filled with dielectric [2^. A cylindrical capacitor was suspended on torsional levers 
inside some external magnetic field parallel to the axis of the cylinder. By applying voltage 
to the capacitor's walls a radial electric field was created. These independent fields possessed 
non-vanishing angular momentum, which was compensated by the motion of the capacitor 
itself. The observed j3| x x pv oc {e — 1) Iq, where Iq corresponds to the e = 1 case, follows 
from (IHdQ- 

Ashkin and Dziedzic observed that the liquid interface bends outwards the liquid in both 

n 

cases when light enters and leaves the liquid[2y]. Contrary to their measurement, the con- 
servation law predicts inward bending. Loudon recently arrived to the same conclusion 
by quantum analysis of the Lorentz force j^. The results of ref. j26i were explained by the 
influence of ponderomotive forces j^], caused by strong focusing of the light in this exper- 
iment. These forces, also used for optical tweezers, are much stronger than contributions 
from the change of the radiation momentum on the boundary. 



5 



Radiation forces can be caused even by redistribution of the energy between quantum 
vacuum and matter. Attraction of two parallel metal plates in vacuum was predicted by 



Casimir 



28| and experimentally observed by Lamoreaux|29|. Electromagnetic field possesses 



finite energy even in the ground state, similar to the quantum harmonic oscillator. The 
presence of dielectric or metallic objects in space alternates the eigenstates of the electro- 
magnetic field. The energy of such system depends on the specific arrangement of the objects 
and some rearrangement can be energetically favorable. Between two metal plates the low 
frequency electromagnetic modes are cut off by the boundary conditions. The smaller the 
separation the smaller is the effective energy of the system, consequently the plates attract 
each other. However, in Casimir case, no momentum is gained by the plate's center of the 
mass according to symmetry considerations. Moreover, to the best of our knowledge, the 
transfer of finite momentum from vacuum modes to matter was not considered yet. 

The zero fluctuations contribution to the equation of motion (fT^ can be expected, since 
the moment of electromagnetic field, similar to its energy, is a quadratic function of E and B. 
Vacuum contribution can not occur, neither in time even media nor in spatially symmetrical 
time-odd (Faraday) materials, due to the self compensation of counterpropagating modes. 
Therefore both time and spatial asymmetries are required. 



The break o: 
materials 



both spatial and time symmetries occur naturally in magnetoelectric 



l3ll |. The dispersion relations for magnetoelectrics are: 
D = eE + xH,B = ilH + x^E. 



(15) 



The same dispersion can be created artificially by applying external electric and magnetic 
fields [32]. In this case, the dielectric properties of the medium e", /i and x depend on the 
external fields E^xt and Bext- For th e sp ecific case of perpendicular electric and magnetic 



fields acting on isotropic material 



For light propagating alon 
equations ((T]) one obtains 



e 
0^0 
e 



'^yu ^ 



X 



Xyx 





(16) 



/i 

y /iy 

Eext X Bext direction, substituting (fT3|) and (fTB|) into Maxwell 
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£^ - Xyx- (17) 



and corresponding modes {E^, Ey, B^, By): 

(1, 0, 0, ^) , (0, 1, -v/^, 0) , (1, 0, 0, , (0, 1, 0) . (18) 



In the case of magnetoelectrics (jl5|) . the term j^Bx^E must be added to the Lagrangian 
(jni). Using ((Zj) one obtains: 

Lme = Lfm + (-Bx^E + —Bx^ x 5) + — (E x t;) e\ . (19) 



Equations of motion (j3|) correspond to the dispersion relations (|T5|l in moving media, while 
fT!?|) becomes: 

p% = ^ ('-^E XB + -EX ifE) - -B X {xB)] • (20) 

The non-compensating moment of a pair of counterpropagating modes in z direction is 
= iXxy — Xyx) (1 + ^jj) I {2n fic). It is obtained by substitution of (fTHj) into (pUj). Taking 
into account all contributing modes and Ap {6) = ApcosO, we obtain: 

TT/2,2TT/uJcut 

i^i+£^ rr ^^^^^Qp-^2dk ^21) 

0,00 

where An = {xxy — Xyx)- The vacuum E^^^ = hLj/2, thus becomes: 

- ' ^-Anl±^^. (22) 
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This expression is significantly different from the Casimir effect, since it is powered by the 
high frequency cut-off. The latter makes it more similar to the Lamb shift phenomenon. 

This effect (I22j) can be evaluated quantitatively by the estimation of the value of An from 
the known experimental data. In crossed external magnetic Bext and electric Ee^t fields. An 
is proportional to magnetoelectric susceptibility 3 1 |3q: 

An^(^P.-lp;)EM^^P.EM (23) 

where Iq ~ 0.3nm is the characteristic interatomic distance. This result follows 
from the spherically symmetric system's fourth order energy terms L = l/AjS^E'^B'^ + 
1/4 (/5|| -3j) {EBy and D = dL/dE relation. The Xxy + Xyx ^ le - 11 was recently 
observed 32| by measurement of magnetoelectric linear birefringence ()17|) in external elec- 
tric Eext = ^ebV/m and magnetic B^xt = 17T fields. The contribution of the non-local 



terms |36l||32j| to An, leading to An oc 1/A, can significantly increase the value of 
However, taking into account that An ^ Xxy — Xyx ~ Xxy +JCyx, the P± ~ O.la.w., which 
corresponds to the experimentally observed An le — llj32| according to (j23|) . is in the 



37| 38|. Therefore An is assumed to be constant in the (j2ip 



range of theoretical predictions 
integral. 

According to ()22|) v^c ~ 50nm/s in external fields Eext = le5V/m and Bext = 17T. 
The cut-off frequency uJcut was chosen to correspond to a wavelength A = 2txc/uj ~ O.lnm, 
since for higher frequencies the molecular polarization vanishes. Density p ~ le3kg/m^, 
An ^ le — 11 and dielectric constant e ~ 1.5 were assumed. The contribution of the 
static field (jl4|) is Vstat ~ 20nm/s. In JRL experiment the estimated velocity from (|14p 
was viaser ~ Ic — 15m/ s (the laser beam intensity was about lebW/rn?). Velocity v^ac is 
not only greater than Vstat and via^, but may also possess the opposite sign. The latter 
follows from substitution of ref. j37| results into (j^Hj) . The experimental measurement of 
p2|) requires effective homogeneity of the matter. Otherwise eqs. (|THl and ()20|1 are not 
valid. Thus the region of the crossed fields must be produced locally, similar to the laser 
beam in JRL experiment. It can be done for instance by immersing a capacitor's electrodes 
inside the liquid. 

Any liquid possess viscosity, causing decay in time of any artificially created flow. The 
dissipation processes can be taken into account empirically by adding correspondent viscous 
forces TjAv to the equation of motion 



p— = r]Av (24) 

while (I14|l defines the boundary conditions. According to |T4jl . the Gaussian laser beam in 
JRL experiment generated Gaussian flow. Eq. ()24|1 . in such case, becomes an ordinary diffu- 
sion equation. The initially gained momentum of the liquid "diffuses" in a plane perpendicu- 
lar to the propagation direction. Thus the characteristic decay time is td ~ ImirrorP/v ~ 
where the mirror's size Imirror ~ 0.25cm and the density to viscosity ratio p/r] 2{){]s / crn^^. 
Radiation pressure in liquids was studied on the time scales comparable with td|8]. Longer 
timescale observations can separate the electromagnetic and liquid flow contributions to the 
measured force, because the measured momentum (fTT|l . in the case of static liquid, corre- 
sponds to Abraham's momentum. It is also feasible in an experiments with strong overlap 
between incident and reflected beams, e.g. measurement in a tube. The motion of matter 
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can not occur under symmetric illumination, thus the measured force must be independent 
of the liquid's refractive index. 

In conclusion relativistic formalism was applied for light-matter Lagrangian derivation. 
Equations of motion were obtained and their correspondence to the Abraham-Minkowski 
controversy related experimental data was demonstrated. The received results correspond 
to Abraham's predictions, while Minkowski's momentum can be obtained from Q without 
its last "motion of the matter" term. Therefore the origin of the controversy lies in the 
underestimation of the fact that the field-matter interaction is impossible without the motion 
of the latter. The vacuum fluctuations induced flow in dielectric hquids with v^ac ~ 50nm/s 
was predicted in external crossed electric and magnetic fields. The significant property of 
this phenomenon is the high frequency vacuum modes contribution, similar to the Lamb shift 
effect. It can be used in future as an investigating tool for zero fluctuations. Other possible 
applications lie in fields of microfiuidics or precise positioning of microobjects, e.g. cold 
atoms or molecules. Initial experimental verifications can be based on artificially created 
random fields. 

The helpful discussions with A. Englander, B. Sfez and Y. Siberberg are gratefully ac- 
knowledged. 
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